Dynamical kickback and noncommuting impurities in a spin chain 
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Abstract 

In an interacting continuous time quantum walk, while the walker (the cursor) is moving on 
a graph, computational primitives (unitary operators associated with the edges) are applied 
to ancillary qubits (the register). The model with one walker was originally proposed by 
R. Feynman, who thus anticipated many features of the Continuous Time Quantum Walk 
(CTWQ) computing paradigm. In this note we examine the behaviour of an interacting 
CTQW with two walkers and examine the interaction of the walkers with noncommuting 
primitives. We endow such a walk with a notion of trajectory, in the sense of sample path of 
an associated Markov process, in order to use such notions as sojourn time and first passage 
time as heuristic tools for gaining intuition about its behaviour. 

Keywords: Continuous time quantum walks; birth and death processes; quantum annealing. 



1 Introduction 

We consider a collection of spin 1/2 systems r(j) = (ti(j), T2(j), Ts(j)), j € A s = {l,2,...,s}, 
coupled with an ancilla qubit a = {oi, 02, 03} through a Hamiltonian of the form: 

1 8-1 

H(a,b) = --^U x ®t + {x + 1) t^(x) + U- 1 ®T + (l) t_(s + 1), (1) 
1=1 

where T±(j) = (ti(j') ± u"2(j))/2. The integers a,b are supposed to satisfy 1 < a < b < s; the 
unitary operators U x act on the state space of the ancilla qubit. We will take, in this note, 

U a = a u U b = a 3 (2) 

and will suppose that all the remaining U x are the identity operator. We will consider an initial 
condition in the eigenspace belonging to the eigenvalue 2 of the conserved number operator 

^^i+^w (3) 

x=l 

and we will refer the system to the orthonormal basis | (xi,X2),( ), where ( G { — 1,1} and 
1 < xi < X2 < s, formed by the simultaneous eigenstates of 73 (x), x £ A s , and 03, belonging, 



respectively, to the eigenvalue +1 of t^{x\) and Tz(x 2 ), to the eigenvalue —1 of the remaining 
t%[x) and to the eigenvalue +1 of Q. 

We will look at the above system from two points of view: 

i. as an Anderson model [TJ with noise on the hopping parameters relative to the links {a, o+l} 
and {b,b + 1}, with the peculiarity that the "random values" of these parameters are 
determined by the non commuting observables o~\ and 03; 

ii. as a version of Feynman's model of a quantum computer j2], where the motion of spin-up 
excitations of the r field (the clock) administers the primitives U a and Ub to the ancilla (the 
register). 

For an extensive analysis of related models in the subspace belonging to the eigenvalue 1 of N3 
we refer the reader to [3]. In the N3 = 1 subspace (because of Peres' conservation laws [I]) the 
presence of an ancilla cannot affect the motion of the clocking excitation. In the N3 = 2 subspace 
we will, on the contrary, give evidence of a peculiar three-body effect involving the two clocking 
excitations and the ancilla qubit, related to the fact that U a and Ub do not commute. 

2 Dynamical kickback 

The study of the system introduced in the previous section is made easy by the fact that the two 
following projectors are constants of motion: 

P± = J2 I (xi,x 2 ),±(-l) e{xi - a)+d{x2 ' a) ) 

( (x 1 ,X2),±{-l) e(xi - a)+9{x2 - a) I, (4) 

where 6{x) = if x > then 1, else 0. If, as we will always do in this note, we consider the 
evolution of the system from the initial state 

|Vo) = 1(1,2), +1), (5) 

we will be interested only in the matrix elements of the Hamiltonian ff (a, b) between states 
belonging to the range of P+: 

h + ((x 1 ,x 2 ), (yi,y 2 )) = 

= ( (x u x 2 ), (_l)«(*i-«O+0(*3-«) |ff (o, 6)| ( yi ,y 2 ), ) = 

= — ^ \^xi,yi {$x2,y2+l + 3x2,y2~l) + ^2, 2/2(^1, yi+1 + ^xi,yi— l)) + 

+ (1 - 9(xi - a))5 Xuyi (5 X2yb Syzjt+i + S X2 ,b+i &yi,b)- (6) 

Looking at the model from the point of view (|n|, it performs quite a trivial computational task: 
starting with the two clocking excitations in positions 1 and 2 and with the register "up" , it re- 
turns, if the two clocking excitations are found in the terminal positions s — 1 and s, the register 
again "up". 

In the case b = a + 1 the only track of the fact that, in applying the identity to the register, the 
machine has temporarily flipped it (by applying <ti), can be seen by comparing the corresponding 
probability amplitude with the one for the free case in which all the U x are the identity operator: 
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(a) 



(b) 



Figure 1: Frame (a): s = 7, a = 4, 6 = 5. Frame (b): s = 7, a = 3, b = 5. Solid lines: the 
probability amplitude —ipt{{s — 1, s), 1) as a function of time, under the initial condition ([5]). The 
dashed lines refer to the free case. 




Figure 2: Graphical representation of the Hamiltonian ([6]) on the weighted graph having the 
set of vertices {(x\,X2) £ A s x A s ,l < < X2 < s} with edges between nearest neigh- 
bour sites. Edges to which (|6]) attributes a positive weight are represented by thick lines. 
Frame (a): s = 7, a = 4, 6 = a + 1 = 5. Frame (b): s = 7, a = 3, 6 = a + 2 = 5. 



the inversion of phase ("dynamical kickback") shown in figure figure 1(a) is easily understood by 
thinking that, while the two clocking excitations move to the right, they restore the register into 
the "up" state by applying minus the identity operator to the register in the successive steps 
(from right to left) a^a\a^a\. 

The above simplistic description of the evolution of the "computation" with two clocking exci- 
tations holds only in the particular case b = a + 1 considered up to now, as made clear by the 



example with b = a + 2 shown in figure 1(b) 

By direct inspection of the Hamiltonian (|6j) and, in particular, of the weighted graph on which 
our quantum walk takes place (figure 2), the role of the positions a and b of the two noncom- 
muting impurities U a = u\ and Ub = 03 is easily understood in the context of an interference 
phenomenon. The term 

(1 — 9{x\ — a)) S Xl)Vl (S xij b ^2,6+1 + ^2,6+1 Sy 2 ,b) 

in h + ((x\, X2), (yi, 2/2)) shows that in the situation b = a + 1 of figure [1(a)] all the amplitudes 
ipt(xi, b + 1), 1 < x± < b are phase inverted with respect to the free case evolving according to 
the finite difference Laplacian 

^1,2/1 (^2,2/2+1 + ^2, J/2-l) + ^2,2/2(^1,3/1+1 + ^i,2/i-l))- 
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In the situation 6 = a+2 of figure 1(b) on the contrary, only the amplitudes ^(xi, 6+1), 1 < x\ < 
a are phase inverted and deviations with respect to the free case take place because of interference 
with the uninverted signal ipt{a+l, 6+1): that this interference can be destructive, suppressing the 
probability of the two excitations ever going beyond the noncommuting impurities, is shown by 
comparison of the two frames of figure [I] A notational remark: for a wave function t/)t((xx, x 2 ), () 
in the range of P+ we are suppressing explicit indication of the argument £. 



3 Sample paths 



It is fairly intuitive to attribute the effect shown in figure |T(b) to the fact that, in the case 6 = a+2, 



not only the computational path <73cri<73<7i = —I is available, but also the path a^a^aiai = I, 
corresponding to the fact that the rightmost cursor can wait in 6 for the leftmost cursor to jump 
in a + 1 ( <T\<T\ = / being thus applied to the register) and then both of them can jump to the 
right of 6 (<73<73 = / being thus applied). 

The above intuition (involving correlations between positions of the cursors at different times) 
can be made more precise in terms of the stochastic process (qi(t), q 2 (t)) associated, according to 
the prescription of [5] (as specialized to the present context in [6]) to the time evolution ipt(xi, x 2 ), 
in H+ = range(P+), of the initial condition The transition probability per unit time from 
site (x\,X2) to site (2/1,2/2) is given by 



v t {yi,y 2 \xx,x 2 ) = \h + ((x 1 ,x 2 ),(y 1 ,y 2 ))\ 

Wt{xi,x 2 ) 

• [1 +sm(Arg(ip t (xi,x 2 )) - Arg(ip t (yi,y 2 )) + Arg(h + ((xi,x 2 ),(yi,y 2 ))))]. 



^(2/1,2/2) 



(7) 



We will indicate by (qi(t), q 2 (t)), v®(yi, y 2 \xi, x 2 ), ■ ■ ■ , the analogously defined quantities in the 
absence of interaction. Simulation of processes with the above transition fields are performed, in 
what follows, according to the first order algorithm outlined in [6]. 

The free process (gi(t), q 2 (t)) (some sample paths of which are shown in the insets of figure |3| 
can, in the region {(xi,x 2 ) £ A s x A s : 1 < x% < x 2 < s}, be described in much the same way 
as the paradigmatic example of |6]: each component starts as a pure birth process (only steps 
to the right {x\ — ► x\ + 1) or upwards (x 2 — > x 2 + 1) are allowed in an initial time interval); at 
each instant each link (edge between nearest neighbour lattice sites) can be traversed only in one 
direction; the allowed direction along a link is inverted each time the probability mass at one 
vertex of the link vanishes. 

Mastering the slalom at (a + 1,6) in the interacting shown by figure |3j requires subtle 

time correlations. There is an instant of time at which most of the trajectories that hit (a + 1, 6) 
are there simultaneously. At a later time they radiate from (a + 1, 6) in many different directions. 
As figure [4(b)] shows, given that the process hits (a + 1,6), it stays there, on the average, for a 
longer time in the interacting case than in the free case. 



4 Conclusions and outlook 

We have tried, in this note, to contribute to the effort of looking at quantum mechanics as a 
source of metaphors suggesting Markov processes with interesting dynamical behaviour, interest- 
ing from the point of view of, say, efficiently crossing (in the sense of [7] ) a graph or a decision 
tree, or sampling a function to be minimized [H |9l [TO] . 

In [6] we have, in this spirit, shown how to mimic, by a Markov process of the class proposed in 
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Figure 3: s = 25, a = 11, 6 = a + 2 = 13. Frames (a) and (b): sample paths of the interacting 
processes qi(t) and q% (t) respectively. Only the trajectories, in our sample of size 10 4 , that hit 
the site (a + 1,6) are shown. In the insets, the corresponding trajectories of the free process 
(q®(t), q^it)) are shown for comparison purposes. 



[5], the diffraction effect due to a sharp initial position. Here we have tried to formulate in the 
same stochastic language the interference effects due to different localizations, as shown in figure 
[2] of discontinuities of the Hamiltonian. 

In the process of doing so we have explored the notions of first passage and sojourn times for 
a quantum walk, which might prove useful from the point of view of suggesting heuristics of 
quantum algorithms, in a context, such as Feynman's [2], in which timing and synchronization 
issues play a major role. 

From the point of view of physics, our analysis raises the question of finding, in the quantum me- 
chanics of Anderson localization by non commuting impurities, an analog of the time dependent 
phenomenon shown by the stochastic process in figures [3] and |4j the sudden formation (see figure 



4(a)), in the situation of figure 2(b) of a probability bubble at (a + 1,6) and its delayed (see 



figure [4(b)] ) bursting in random directions. 
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(a) 



(b) 



Figure 4: s = 25, a = 11, b = a + 2 = 13. Conditional cumulative distribution function of the 
first passage time (Frame (a)) and sojourn time (Frame (b)) at site (a + 1,6), given that in the 
time interval (0, s) the process visits (a + 1, b). Solid line: the interacting case. Dashed line: the 
free case. Both distribution functions are estimated from the subsample of those trajectories that 
do hit (a + l,b). 
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